Electronic properties of EuBg in the ferromagnetic regime: Half-metal versus 

semiconductor 
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To understand the halfmetallic ferromagnet EuBe we use the Kondo lattice model for valence and 
conduction band. By means of a recently developed many-body theory we calculate the electronic 
properties in the ferromagnetic regime up to the Curie temperature. The decreasing magnetic order 
induces a transition from halfmetallic to semiconducting behavior along with a band broadening. 
We show the temperature dependence of the quasiparticle density of states and the quasiparticle 
dispersion as well as the effective mass, the number of carriers and the plasma frequency which are 
in good agreement with the experimental data. 

PACS numbers: 71.10.-w, 75.50.Pp, 71.30.+h, 75.10.-b, 75.3G.-m 



I. INTRODUCTION 

In the last decade EuBg with its interesting and still 
not fully understood magnetic and electronic properties 
has attracted considerable interest. For low temperatures 
it is a ferromagnetic semimetal becoming a paramagnetic 
semiconductor for higher temperatures. It crystallizes 
in the cubic CaBg structure (spacegroup 221, Pm3m). 
The divalent Europium ions occupy the corners of the 
cubic unit cell and at the body centered position sits an 
octahedron consisting of six Boron atoms. The lattice 
constant a=4.19A does not change with temperatureii 

EuBg exhibits two consecutive phase transitions at 
Tc=12.5K and Ta/ = 15.5K:1 T^ is the bulk Curie tem- 
perature. The magnetization follows a Brillouin func- 
tion and the saturation value is 7fiB due to the exactly 
half- filled 4f shell in Europium (2i Nonetheless there are 
evidences of phase separation above Tc up to SOK^, 
where bound magnetic polarong- are responsible for the 
cusplike peak in resistivity and the colossal magnetore- 
sistance at Tm<^ Concomitant with the bulk magnetic 
phase transition at Tc is a transition from semimetal to 
semiconductor . ^' ^ 

For low temperatures the semimetallic behavior is 
proven theoretically and experimentally. Early local den- 
sity approximation (LDA)-- as well as recent band struc- 
ture calculations in LDA-|-U^ show an overlap of valence 
and conduction band around the X-point of the Brillouin 
zone. Shubnikov - de Haas and de Haas - van Alphen 
measurementsJJ^jii agree with that picture as they report 
electron and hole pockets on the Fermi surface at the 
X-point. The carrier density is about 0.01 per unit cell 
for low and 0.001 for higher temperatures:^^ Due to the 
known lack of correct predictions of band gaps in the 
LDA the authors in 9] suggest EuBg to be halfmetallic 
in the ferromagnetically ordered phase. That scenario is 
not ruled out by experiments albeit not proven also. 

The resistivity shows a strong peak at Tm and also 
indicates a metal - insulator transition. Up to Tc and 
Tm an increase of the carrier concentration, a reduction 
of the effective mass or a combination of both is respon- 
sible for a strong red shift of the plasma frequencyii^ The 



specific electronic properties seem to stem from tempera- 
ture dependent changes of valence and conduction bands 
induced by the ferromagnetic alignment of the local mo- 
ments. The other way around the ferromagnetism is at- 
tributed to an RKKY mechanism where the itinerant car- 
riers are responsible for the magnetic exchange. It looks 
as if the system will be best described by the two band 
Kondo lattice model introduced in the next section. 

Similar models but with focus on other properties of 
EuBg are already studied in a mean-field treatmenti^ 
and in dynamical mean-field theory. The spin polaron 
modeli^ explains the behavior in resistivity and the colos- 
sal magnetoresistance effect around T m ■ This paper will 
focus on the range from zero temperature to Tc where 
the magnetic and metallic phase transitions take place. 



II. THEORY 

For modeling EuBg we will use the Kondo-lattice 
modeljii also known as s-f model. The spins of the lo- 
calized 4f electrons result in a local magnetic moment 
which interacts with the spin of the itinerant electrons in 
the valence(v) and conduction(c) band, respectively. The 
complete two-band Hamiltonian in second quantization 
is 
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The first term in describes the kinetic part of the elec- 
trons in band b with the hopping integrals T^^ between 
two lattice sites and that are connected to the 
wave vector k dependent Bloch energies via Fourier 
transformation T^*^ = e^e*''^'^ . The chemical po- 

tential is ^ and aJjo.(afcjcr) the creation (annihilation) op- 
erator of an electron with spin a at lattice site Rj in 
band b. 
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The second term in describes the s-f interaction be- 
tween the itinerant band electrons and the local moments 
Sj. Taai are the Pauli matrices and j'"'" the coupling 
constants. In contrast to the kinetic part the interaction 
term does not need to be diagonal in band indices. It 
might be possible that during the s-f interaction the itin- 
erant electron in band b hops to the other band b'. That 
possibility is included in (Q. 

In EuBg the valence band is formed from Boron 2p 
and Europium 4f states and the conduction band from 
Boron 2p and Europium 5d states which overlap at the X- 
point^. The possible interband transitions are negligible 
by symmetry reasons. 

Therefore we will set j'''" diagonal, resulting in the 
following Hamiltonian 



For some special cases the equation of motion can be 
solved exactly but for most of the cases one needs to 
apply specific approximations for the higher Green func- 
tions to decouple. Decoupling means it is possible to 
introduce a self energy Y]\_^{E) which formally solves 



and thus the Green function reads 
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B. Self Energy 
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Iri is the kinetic term for band b 
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describes the s-f interaction for band b 
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where we used ribja- = a\,j„o,bjo 



Scr] - 5cri and S". 



Sj + iZaSj . The first term in Q describes an Ising-like 
interaction and the second the spin-flip processes. 

The coupling constant J** is positive (ferromagnetic 
alignment) for the conduction band electrons and nega- 
tive (anti-ferromagnetic alignment) for the valence band 
electrons^. We assume the same magnitude J for both 
bands, hence J"=-J and J'^=J. 

It goes without saying, that the k-independent cou- 
pling constant is an oversimplification. However, this is 
unavoidable due to mathematical difficulties. 



A. Green Function 

We are mainly interested in the single electron Green 
function 



GijE) = ((a,k.;aL)> 
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from which all electronic properties can be calculated. 
The operator ajj.^(abkcr) is the creation (annihilation) op- 
erator in k-space. To derive the Green function we solve 
the equation of motion 

EGi,{E) = h+{{[a,k,,H'^ + H'A;alj. (7) 



In this section we will introduce an appropriate ap- 
proximation of the self energy for the almost empty con- 
duction band, which is ferromagnetically coupled to the 
local magnetic moments, and the almost filled valence 
band which is anti-ferromagnetically coupled. 

For the conduction band it is possible to apply the 
interpolation formula for the self-energy given in [2(il | 
which is valid for low densities (e^O) and covers the fol- 
lowing exactly solvable special cases: 

• Zero-bandwidth limit (sk Tq) 

• Ferromagnetically saturated semiconductor 
i{S^)=S, e=0) 

• Second order perturbation theory in J 

• High energy expansion 
The resulting self-energy is 
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where = 



= S{S -I- 1) — m„{mcr -t- 1) and 
k E-e-k+tJ. ^"^^^ propagator. The mag- 



netization {S^) will be calculated via the Brillouin func- 
tion which is in agreement with the measurements in 0] . 

Since the valence band is almost completely filled, let 
us change for one moment from the electronic picture to 
a hole like quasiparticle picture. Such a hole corresponds 
to a missing electron in an otherwise completely filled va- 
lence band. To keep the negative coupling constant for 
the electrons in the valence band, the coupling constant 
for the holes must have opposite sign (J^(e) —J'"{h)). 
The self energy for the holes in the low density approxi- 
mation (h^O) from '2^ then reads 
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Returning to the electron picture through reversing the 
spin {Mcr{h) —M„{e)), the self energy for the electrons 
in the almost completely filled valence band then is 
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In fact this self-energy covers the above mentioned special 
cases for full bands. 

Now we can put the selfenergies for both bands, H1U|) 
and H12|l in and calculate the Greenfunction. 



III. RESULTS 

As described in the previous section we are able to 
calculate the one electron Green function with the given 
self energies for both bands. In this section our results are 
presented. The band index will be skipped as it should be 
clear that the following calculations are formal equivalent 
for both bands. 

Several free parameters need to be adapted. The spin 
is S=| due to the exactly half filled 4f shell of Europium 
and the coupling constant J=0.2cV is in accordance with 
current literature. We assume a modified tight binding 
dispersion £k for simple cubic^^ with the band gap lo- 
cated at the X-point instead of T. 

£k = To \ cos{k.xa — tt) + cos{kya) + cos{kza)j (13) 

Here Tq is the band's center of gravity and W the band- 
width. The widths of the tight binding bands are cru- 
cial for our model calculation since the effective coupling 
strength is The mean values of the anisotropic ef- 
fective masses m* given in Q are 1.2me and 0.35me for 
valence and conduction band, respectively. From that we 
can calculate the corresponding bandwidths 



W = 



(14) 



which leads to W"— 2.2cV and W=7.4cV and is in agree- 
ment with the approximated widths from the band struc- 
ture calculations for the hexaboridsSiS. In accordance 
with the measured band gap for the non-magnetic (J=0) 
CaBg in |2ll we assume an intrinsic gap of A=:0.25eV. 
The chemical potential is determined to result in compen- 
sated electrons and holes since the intrinsic band struc- 
ture is semiconducting. 



A. Q-DOS 

Figure n shows the quasiparticle density of states (Q- 
DOS) for the electrons in the valence and conduction 
band at different temperatures (OK, |Tc, Tc) 
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FIG. 1: Quasiparticle density of states for spin down (left) 
and up (right) electrons in valence (straight line) and conduc- 
tion band (dashed line) at T=OK (top), |Tc (middle) and 
To (bottom), S=3.5eV, J=0.2eV, W"=2.2eV, W'^=7.4eV, 
A=0.25eV. 



Since JS^W we are in the weak coupling regime, al- 
though the effective coupling strength ^ is different for 
conduction and valence band which leads to different re- 
sults. 

For the conduction band the Q-DOS is shifted with 
magnetization in a mean-field like manner by ^JM„. 
Even so one can see an overall band broadening with 
decreasing magnetization. 

For the valence band some correlation effects are vis- 
ible. At T=OK the spin down band is rigidly shifted by 
iJS. Since the local moments arc completely aligned no 
spin exchange is possible for spin down electrons (keep 
the anti-ferromagnetic coupling in mind) hence there is 
no deformation of that band. In the spin up channel one 
sees a dip between two overlapping subbands. The lower 
subband describes electrons which flip their spin during 
the interaction with the local spin system, emitting a 
magnon. Thus their lifetime is decreased. The upper 
subband is connected to a stable quasiparticle known as 
magnetic polaron which propagates through the lattice 
dressed by a virtual cloud of magnons (for more details 
see [23). 

At higher temperatures the scattering for both spin 
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FIG. 2: Quasiparticle spectral density at the F-point for spin 
down quasiparticles in the valence band at T=OK, \Tc and 
Tc, S=3.5eV, J=0.2eV, =2.2eV, =7.4eV, A=0.25eV 



orientations becomes more likely, explaining the defor- 
mation of the Q-DOS in both spin channels. With de- 
creasing magnetization (increasing temperature) the dif- 
ference between the up and down spectrum decreases and 
vanishes at Tc due to the missing alignment of the local 
moments. 

What is also shown in figure ^ is a transition from 
halfmetallic to semiconducting behavior undergoing the 
magnetic phase transition. The overlap of the spin up 
bands at OK is 0.5eV and the gap between valence and 
conduction band at Tc is 0.04eV. The band gap at Tc 
is in accordance with tunnel experiments^. 

With a mean-field treatment of the Hamiltonian one 
could also show the halfmetal to semiconductor tran- 
sition. The tight binding dispersion would be shifted 
rigidly by ^JM„, resulting in an overlap at T==OK of 
0.45eV for spin up and the band gap at Tc of 0.25eV. 
The smaller band gap within our approach results from 
the interaction between the electrons and the local mo- 
ments even above Tc. Although the net magnetization 
at Tc is zero, the statistically distributed magnetic mo- 
ments influence the electronic structure which results in 
the band broadening. In a meanfield treatment the in- 
teraction is simply switched off at Tc- 



B. Spectral Density 



The spectral density 
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can be compared directly with angle and spin resolved 
(inverse) photo emission experiments. 

Figure 13 shows the spectral density at r=(0,0,0) for a 
spin down quasiparticle in the valence band at different 
temperatures. As explained in the previous section, at 
zero temperature the spin down electron in the valence 
band cannot exchange its spin with the local moments 
because they are fully aligned. Thus the quasiparticle 
has infinite lifetime. That fact is visualized in figure 
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FIG. 3: Spectral density as grayscale at k-points along the 
irreducible Brillouin Zone for spin down (left) and up (right) 
quasiparticles at OK (top), |Tc (middle) and Tc (bottom), 
S=3.5eV, J=0.2eV, W =2.2eV, =7.4eV, A=0.25eV 



El as a sharp peak for T=OK. For higher temperatures 
and less magnetic ordering, magnons are excited and can 
be absorbed by the spin down electron in the valence 
band. Therefore the lifetime is decreased with temper- 
ature which is represented by the wide peaks in figure 

m 

That evaluation of the spectral density could be done 
for every k-point. In figure|3|we show the spectral density 
as grayscale along the standard symmetry points in the 
Brillouin Zone, which makes a temperature dependent 
band structure. 

The same effects we noticed in the Q-DOS (Fig. ^ can 
be seen in the density plots (Fig. EJl. The spin up con- 
duction band and spin down valence band is undeformed 
and rigidly shifted, represented by the thin lines. The 
smearing occurs where spin flip scattering is possible. 

For example in the plot for spin up quasiparticles in 
the valence band at T=OK (fig. 13 top right, lower band) 
we can distinguish between an undeformed, sharp region 
around the X-point and a blurred region. These are con- 
nected to the two subbands mentioned in the previous 
section. The lower one describes quasiparticles that are 
involved in spin fiip scattering, their lifetime is decreased. 
Whereas for quasiparticles in the upper subband no real 
spin flip processes are possible because there are no states 
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FIG. 4: Temperature dependence of the effective mass m 
(dashed Une) and occupation number n (straight line) of spin 
up quasiparticles in the conduction band, S=3.5eV, J=0.2eV, 
=2.2eV, W= =7.4eV, A=0.25eV. 




FIG. 5: Temperature dependence of the plasma frequency for 
the conduction band. 



in the corresponding spin down region. Their hfetime is 
infinite. 

At higher temperatures the spin flip processes are pos- 
sible for spin up and down. It occurs especially in the flat 
regions of the dispersion where the possibility of scatter- 
ing is enhanced due to the smaller velocity. 

Along the T — R direction occurs a band splitting ap- 
proaching Tc- In a photo emission experiment it would 
be difficult to resolve the light shaded regions in figure 13 
because of the short lifetime of the quasiparticles. There- 
fore one might wonder about two bands coming from F 
and R which do not fit together. They can be explained 
with our many-body theory of the interaction between 
the electrons and the local moments. 

The effects in the conduction band are less pronounced 
due to the weaker effective coupling. Nevertheless one 
can see that the lines get smeared with temperature due 
to the interaction. 



C. Plasma Frequency 

Now we will show the qualitative behavior of the ef- 
fective mass, occupation number and plasma frequency. 
That can be compared directly with the experiments and 
will be a measure whether our Kondo lattice model is able 
to explain some of the properties of EuBg. 

Via the partial derivative of the real part of the self 
energy R^- with respect to the resonance energies ii^ko- — 
£k — /i + ^CT(-EkCT) it is possible to derive the effective mass 
of the quasi-particles 



(17) 



The effective mass at the X-point for a spin up quasipar- 
ticle in the conduction band is shown in figure 0] (dashed 
line). One can see a slight increase with increasing tem- 
perature. That is in qualitative agreement with the ex- 
perimental results in |0| and fl^ . 

The occupation number can be computed by 



nAT) = E(«L«k.)(T) (18) 

k 

= J f^{T,E)gAT,E)dE (19) 

where /_ (T, E) = (l+exp{^^^)) ^ is the Fermi function 
(ks is the Boltzmann constant). 

Figure ^ (straight line) shows the occupation number 
for spin up quasi-particles versus temperature. In the 
previous section the decreasing band overlap with tem- 
perature due to the reduced magnetization was intro- 
duced. This is why the occupation number decreases 
dramatically by an order of magnitude. This agrees with 
the Hall data in [Tl |. 

FigureElshows the calculated plasma frequency versus 
temperature 



eom(T) 



(20) 



with the electric charge e and the dielectric constant eg. 
The reduction of free carriers and the increasing effective 
mass with increasing temperature lead to a red shift of 
the plasma frequency. That is in excellent agreement 
with magneto optical results in [l^. 



D. Half metalicity 

As already mentioned we derive half metallic behavior 
in the whole ferromagnetic regime. To clarify that, figure 
|H1 shows the shrinking band overlap for spin up (straight 
line) and the gap for spin down quasiparticles (dashed 
line) versus temperature. At zero temperature the spin 
up bands overlap by 0.5eV and the spin down bands are 
separated by 0.9eV. With decreasing magnetization the 
spin up and down bands are shifted and merge at the 
Curie temperature, resulting in a band gap of 0.04eV at 
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FIG. 6: Relative temperature dependence of the band overlap 
(straight line) for spin up and the band gap (dashed line) for 
spin down quasiparticles. 

IV. CONCLUSIONS 

Starting from the two band Kondo lattice Hamiltonian 
with opposite sign of the coupUng constants for valence 
and conduction band electrons we were able to reproduce 
the metal-insulator transition for EuBg concomitant with 
the magnetic phase transition. The synchronous band 
broadening is verified through our model, too. The ex- 
perimentally measured dependence of effective mass, oc- 



cupation number and plasma frequency are qualitatively 
confirmed. 

In our model the change in the carrier concentration 
dominates over the change of the effective mass, whereas 
in the model proposed by Hirsch^'' only the change of the 
effective mass is made responsible for the specific prop- 
erties of EuBg. 

With the chosen parameters, motivated by exper- 
iments, EuBg turns out to be not only metallic 
but halfmetallic in the ferromagnetic regime. That 
would need further experimental investigations since the 
halfmetallic behavior is announced as a possible result of 
future GW-calculations in 9] and not explicitly excluded 
in the Shubnikov - de Haas and de Haas - van Alphen 
measurements. 

A possible refinement of our calculation will be to con- 
sider a real band structure for EuBg. Above Tc up to 
30K is a regime of phase separation where bound mag- 
netic polarons and percolation effects play a significant 
role. That is not included in our model yet, since the 
magnetization is not calculated self consistently. That 
should be done in future corresponding to j22| |. Further 
on it would be interesting to study the diluted system 
Caa;Eui_j;Bg as was done in p^ . 
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